Majorana fermions in superconducting nanowires without spin-orbit coupling 
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We study nanowires with proximity-induced s-wave superconducting pairing in an external mag- 
netic field that rotates along the wire. Such a system is equivalent to nanowires with Rashba-type 
spin-orbit coupling, with strength proportional to the derivative of the field angle. For realistic 
parameters, we demonstrate that a set of permanent magnets can bring a nearby nanowire into the 
topologically non-trivial phase with localized Majorana modes at its ends. This occurs even for a 
magnetic field configuration with nodes along the wire and alternating sign of the effective Rashba 
coupling. 



Currently, there is an intensive search for materials and 
geometries that can support topological states, such as 
Majorana bound states appearing in systems with a p- 
wave superconducting order. One motivation for this ef- 
fort is the possible applications of Majorana states as the 
basis for topological quantum computing,^ for example in 
hybrid structures combined with other qubit systems^^ 
to achieve universal quantum computing. This large ac- 
tivity is inspired by the recent suggestions to use an ordi- 
nary s-wave superconductor in proximity to a topological 
insulator^ or to systems with strong spin-orbit coupling 
and a strong Zeeman or exchange field to en gineer super- 
conductors with p-wave order parameter pHlSUMUl The 
role of the spin-orbit interaction is to mix the spin direc- 
tions, polarized by the applied field, such that there is an 
overlap with the induced s-wave superconducting order. 
For example, a one-dimensional system with spin-orbit 
coupling and a Zeeman field perpendicular to the spin- 
orbit field in proximity to a s-wave superconductor^ ef- 
fectively (when projected to the lowest band) reduces to 
Kitaev's model™ known to support Majorana modes at 
its ends. The hitherto proposed systems, which combine 
s-wave superconductors, Rashba spin-orbit coupling and 
splitting of spin degeneracy, include heterostructures of 
superconductors and ferromagnets, with spin -orbit cou- 
pling at the superconductor surface ^ 17 * 19 * 21 * or an addi- 
tional layer of strong spin-orbit coupling s emicond uctors 
with proximity induced superconductivity ! 15 * 16 * 1 ^ A one- 
dimensional version of this allows the ferromagnet to be 
replaced by an external magnetic field IHUJ] Other recent 
one-dimensional suggestions included using nanotube^l 
or chains of quantum dotsP^* 

The realization of Majorana fermions could be easier 
if one could relax on the requirement of spin-orbit inter- 
action in the material. Braunecker et aZpS showed that 
in one dimension a spiralling magnetic field is equivalent 
to a Rashba spin-orbit coupling in the local frame. Such 
a situation can be realized by a helical nuclear magnetic 
orderP^ Along the same lines, Choy et al^ proposed to 
use an array of magnetic nanoparticles to create a sim- 
ilar situation, but in this case the spatial dependence 
of the magnetization direction was due to non-collincar 
alignment of the magnetic moments of the nanoparticles. 




FIG. 1. (Color online) Sketch of a quantum wire subjected 
to a spatially varying magnetic field due to a set of magnetic 
gates deposited on the substrate. 

Similarly, it was argued that a superconductor with he- 
lical magnetisirPH produces a p-wave superconductor. 

Here, we investigate a different approach based on spa- 
tially varying magnetic fields generated by easily fabri- 
cated magnets of submicrometer sizes. To demonstrate 
the principle, we focus in our calculations on the geome- 
try shown in Fig. 1, namely, a one-dimensional wire next 
to a set of parallel permanent magnets, with either alter- 
nating or parallel magnetizations [see also Fig. 2(a)]. Al- 
ternatively, the spatial variations in the Zeeman field di- 
rection can be created in a bend wire with an anisotropic 
g-factor and a constant field. 

To model such a geometry, let us consider a one dimen- 
sional system with a magnetic field that changes direction 
along the wire. The wire is placed in tunnel contact with 
a s-wave superconductor, which induces a pairing inter- 
action potential A (assumed real and local below). Using 
the Nambu four-vector basis * = (\I/ t , vf/J, -*|), the 
Hamiltonian for the nanowire is H = § J" d^H^, where 

n = n Q + n s , (i) 

where the normal and superconducting parts are, respec- 
tively, 

Ho = - /.J r 3 + ifl/*flB(0 a, H S = An, (2) 

where B(£) is the effective magnetic field field along the 
wire coordinate £, g is the effective Lande g-factor, Ui 
(i = 1,2,3) are Pauli matrices operating in spin-space, 
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while n are Pauli matrices in electron-hole space. We 
have not included an intrinsic spin-orbit coupling in the 
Hamiltonian ([2]), but it should be noted that in order to 
align with the effective spin-orbit field defined below, an 
intrinsic spin-orbit field should be perpendicular to the 
plane of the applied field. 

Below, we present numerical calculations of the Hamil- 
tonian in Eq. ^ for a magnetic field configuration set up 
by a set of permanent magnets, causing a spiralling di- 
rection of the field. To understand why such a system is 
equivalence to a Rashba-type spin orbit interaction, we 
follow Braunecker et al¥^ and perform a unitary trans- 
formation. We rotate the z-axis of the local spin basis to 
the direction of the magnetic field by the unitary operator 
U — exp(i(ip/2)<r xy ), where a xy — ((B x z) ■ <r)/|B x z| 
and cos ip = B • z/B, with B = |B|. The transformed 
Hamiltonian H = U'HU then becomes 
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while the rotation, of course, leaves the singlet pairing 
interaction invariant Us = Us- The rotation gener- 
ates two new terms, Hr = (H/mi)WU'p^T 3 and %2 = 
{—h 2 /2m)U^U"T 3l where the primes denote differentia- 
tion with respect to £. The term Hr, being proportional 
to the momentum operator, is a Rashba-type spin-orbit 
coupling that reduces to 



Hr = 



dip 



(4) 



This simplifies further if the field lines and the wire lie 
in a single plane, which is indeed the case for the lat- 
eral structures considered here (see Figs. 1 and 2). In 
that situation, we can choose z to be in the plane of the 
field and hence a' xy = 0. To avoid a sign change of a xy 
when B || z, we choose U = exp(iipa±/2), with a± con- 
stant and ip continuous. For the lateral configuration, the 
Rasbha-type spin-orbit coupling part of the Hamiltonian 
thus becomes 



Ur = aoffOjj? 4 T3, 



(5) 



where we define an effective spin-orbit interaction coeffi- 
cient as 



h dip 
2m d£' 



(6) 



The second new term in the transformed Hamiltonian ([3]) 
becomes H 2 = {h 2 /2m)[(ip' /2) 2 - iip"a±/2]r 3 , where the 
first term renormalizes the chemical potential. The last 
term has the form of an imaginary magnetic field parallel 
to the plane of the spin-orbit field, and it appears because 
the total Hamiltonian has to be Hermitian. Alterna- 
tively, this term can be absorbed into Eq. ^ by writing 
it in symmetrized form, i.e., (a c ffP£ + p^a c s)/2. In or- 
der to drive the one-dimensional superconductor into the 
p-wave state, a large spin-orbit coupling and a Zeeman 



* A ; 



~ ¥1 W TTTT " TUT I 



HOIQIIHIOIffl 




2.5 
£ / l^m] 



2.5 

£ / \M m ] 



FIG. 2. (Color online) Sketch of two realistic systems having 
a series of permanent magnets with either alternating (left) 
or parallel (right) magnetization directions. The top panels 
show the field lines and it is clearly seen that the field changes 
direction a number of times along the wire. The bottom pan- 
els show the amplitude of the magnetic field B / Bq in red 
(light gray), as well as the effective induced spin-orbit cou- 
pling Q ff in blue (dark gray) (calculated with effective mass 
m — 0.014m e ). The overall scale for the magnetic field Bq is 
set by the magnetization of the magnets. The magnets have 
widths 600 nm, heights 330 nm, and the gap between them is 
200 nm. For the left configuration the distance to the wire is 
100 nm, and 50 nm for the parallel configuration. 



field perpendicular to the spin-orbit field are desirableP^l 
Therefore the rotation of the magnetic field should be 
optimized to have a large first derivative, <p' . 

An optimal and illustrative example is there- 
fore a sinusoidally rotating magnetic field B(£) = 
i3 c (sin(^/i?)) i 0, cos(£/i?)), which gives the transformed 
HamiltoniarPZl 
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H B B c a 3 . 
(7) 

This model is seen to be identical to that of a one- 
dimensional wire with Rashba spin-orbit coupling, with 
a c ff = h/2mR. To have a large effective spin-orbit cou- 
pling, one should therefore use low mass materials and 
engineer a large curvature for the magnetic field. Using 
parameters relevant for InSb (which also has a large g- 
factor, g « 50), we set m = 0.014m e , which together 
with R = lOOnm gives a c s ~ 3 x 10 4 m/s. This is similar 
to the spin-orbit coupling strength in, for example, InAs. 

A realistic way to generate a magnetic field that ro- 
tates along the wire is to place the wire next to a set 
of permanent magnets, as illustrated in Fig. 1. The 
magnets could, for example, be made of Co. The left 
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configuration in Fig. 2(a) has a staggered configuration 
(which could be realized by the magnets having different 
sizes, so that their hysteresis loops are also different), 
while the configuration to the right has aligned magnets. 
We have computed the fields, with dimensions that are 
easily fabricated with micro technology. For both cases, 
Fig. 2(b) shows the amplitude of the magnetic fields B (in 
units of Bq, set by the permanent magnets) and the ef- 
fective Rashba spin-orbit coupling parameters a c g given 
by Eq. (J6|. The configuration with alternating magne- 
tization directions is seen to give a field close to that 
of the optimal sinusoidal variation (which has constant 
B and a c s). For the configuration with parallel mag- 
netizations, the effective Rashba coupling changes sign 
along the wire, but nevertheless with a sufficient inte- 
grated weight to induce a phase transition to a state with 
Majorana fermions, as shown below. 

We now discuss the conditions for bringing the 
nanowire with proximity induced pairing interaction into 
the topologically non-trivial regime, and thus to have Ma- 
jorana fermions at its ends. Starting with the case of 
a superconducting wire with sinusoidally rotating mag- 
netic field, leading to the normal part in Eq. ([T]), there is 
a transition to a non-trivial phase wherP^ 



gHB\B c \ > VlAP + ^-^/Smi? 2 ) 5 



(8) 



where /x is the chemical potential. It is interesting to 
note that the curvature term shifts the chemical poten- 
tial, allowing for a slightly larger density of electrons in 
the non-trivial phase. 

In the general case, the transition to the non-trivial 
phase cannot be found analytically and we instead resort 
to numerical methods; guided by the above conclusions, 
that it is desirable to have a fast rotation of the direction 
of the magnetic field. In our numerical study we use 
original Hamiltonian in Eqs. and We have studied 
two criteria for the existence of the end modes, namely, 
that (i) a change of the determinant of the zero-energy 
reflection matrix from 1 to — 1 a t the transition to the 
topological non-trivial phase, 24 26 and (ii) the existence of 
a pair of zero modes. We calculate the scattering matrix 
of the wire connected to two single-mode metallic leads 
using the expression^ 



S(Q) 



1 + mvW^U- l W aT 
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where W is the coupling matrix between to the two leads 
(with density of states v) at the end point of a discretized 
version of the Hamiltonian and it is thus a matrix with 
dimension 2 x N, with N being the number of sites in 
the chain, having entries only at the upper left and lower 
right corners. Further, SqSq is the ^-matrix in the ab- 
sence of coupling to the wire. For the numerical calcu- 
lations, the wire is discretized as a tight-binding chain 
with hopping matrix element t between sites separated 
by a and the chemical potential in the bottom of the 
band with approximate quadratic dispersion, so that we 
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FIG. 3. (Color online) Phase diagram for the superconduct- 
ing nanowire subject to the magnetic field in the two configu- 
rations shown in Fig. [2] calculated for 5 — ^m-long wire with 
m = 0.014m e and A = 0.3meV. The topologically non-trivial 
phase occurs when the determinant of the reflection changes 
from 1 to — 1, which is shown in the phase diagram in (a) for 
changing magnetic field B and chemical potential fi and also 
in the lower panels along the dashed (red) line /i = 0.3A. The 
lower panels show the lowest positive eigenvalue Eq, which 
goes to zero after the topological phase transition, as well 
as the next positive eigenvalue E\. After the transition the 
difference between the two is equal to the gap of the topo- 
logically non-trivial state. Both configurations in Fig. 2 are 
seen to exhibit a phase transition. In (b) and (c) the lines 
in red (light gray) show the value of the determinant of the 
reflection matrix. 



can relate to the effective mass via m = f? /2ta. Because 
of finite-size effects the result for the reflection matrix 
determinant is not independent of the choice of is, and 
for the actual calculation we use a value showing most 
clearly the transition. In Fig. 3 we use N = 200 and 
have checked that the results do not change significantly 
by changing N. 

In addition, we have numerically diagonalized the dis- 
cretized version of the Hamiltonian in Eq. ([2| and the 
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FIG. 4. (Color online) The localized Majorana modes and 
their spin texture. The top panel shows the total weight of 
the end Majorana modes, y a and 76, while the lower panel 
depicts the local spin direction for the configuration of alter- 
nating magnetizations (Fig. 2 left). The length of the arrows 
in the inset corresponds to the amplitude of the wave function. 
Notice that the :r-axis of the inset is rescaled. 



topological phase transition is characterized by the point 
where the system supports a zero-energy state, while hav- 
ing a gap for the continuum states. Figure [3] shows the 
results for determinant of the reflection matrix together 
with plots of the lowest eigenenergies and the gap. In 
Fig. |3ja) the determinant is seen to have a sharp transi- 
tion from 1 to — 1. At the same point where this transi- 
tion occurs, a zero mode appears. The zero mode con- 
tinues into the non-trivial phase, while the gap (shown 
as the energy of the next lowest positive eigenenergy) 
increases for larger values. 

Interestingly, the configuration with parallel magnets 
also shows a topological phase transition, with similar 
values of the gap and the necessary magnetic field, see 
Fig. 3(c) (it should be noticed, however, that for this 
case the wire should be closer to the magnets to compen- 
sate for the partial cancellation of the field lines coming 
from neighboring magnets). This is interesting because 
it shows that Majorana fermions can be generated even if 
the magnitude of the magnetic B goes below the critical 
field in a small region of space and even if the effective 
Rashba spin-orbit has an opposite sign in this region-see 
the lower right panel of Fig. 2 for the detailed shapes. 
How robust the topological phase is to having regions 
with subcritical magnetic field is an interesting question 
that deserves further studies. 

Finally, we investigate the structure of the Majorana 
fermions. Each end of the wire supports one Majorana 



state, which we denote j a and %. In the four-spinor 
Nambu basis 7 = (uf, Ui,vi,—Vf), where u and v are the 
electron and hole amplitudes, respectively, a Majorana 
state has — uf and — ulPS However, for our finite 
system the Majorana states are coupled by an overlap ex- 
ponentially small in the distance between them. There- 
fore the numerically determined eigenstates are linear su- 
perpositions of "f a and 7b with a small finite energy. If the 
numerical eigenstates are denoted 71 and 72, we gener- 
ate the states that would correspond to Majorana states 
for an infinitely long wire as 7 a = (71 + 72 VV^ and 
7b = — 72)/v / 2. These states are plotted in Fig. [4] 
The top panel shows the total weight |7 a /h| 2 of the two 
Majorana states localized at the ends of the wire. It is 
also interesting to look at the spin direction of the lo- 
calized Majorana modes, 29 which is also relevant for a 
recent proposal of transfer of quantum information be- 
tween Majorana qubits and spin qubitsP In Fig. [4J we 
plot the local spin direction of the electron part defined 
as (see Ref. \M S a/6 (0 = (l a/b \0(^ <8> |(1 + r z )\ la/h ). 
If one compares the spin direction with the field lines in 
Fig. |2j it can be seen that the direction is dictated by 
the direction of the magnetic field and can therefore be 
tuned. The spin polarization of the Majorana modes re- 
sults in spin-specific tunnel coupling, which can be used 
for both manipulation as in Refs. Eland 0] and detection 
as in Red. El 

In summary, we have shown that quantum wires 
with proximity induced superconductivity and a spatially 
varying magnetic field created by a realistic configura- 
tion of permanent magnets can have topological excita- 
tions in the form of Majorana fermions. This system has 
advantages over proposals that require a specific form 
of Rashba spin-orbit coupling. Moreover, by choosing 
materials with a large (/-factor, a modest field is suffi- 
cient to meet the condition that the Zeeman exceeds the 
induced pairing potential. Finally, we studied the Ma- 
jorana states spin texture, which can be tuned by the 
magnetic field configuration. 
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